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We extend previous analyses of wide-angle correlations in the galaxy power spectrum in redshift 
space to include all general relativistic effects. These general relativistic corrections to the standard 
approach become important on large scales and at high redshifts, and they lead to new terms in the 
wide-angle correlations. We show that in principle the new terms can produce corrections of nearly 
10% on Gpc scales over the usual Newtonian approximation. General relativistic corrections will be 
important for future large-volume surveys such as SKA and Euclid, although the problem of cosmic 
variance will present a challenge in observing this. 



Upcoming surveys of galaxies and HI will probe increasingly large scales, approaching and even exceeding the Hubble 
scale at the survey redshifts. On these cosmological scales, surveys can in principle provide the best constraints on 
dark energy and modified gravity models - and will be able to test general relativity itself. Furthermore, constraints 
on primordial non-Gaussianity from large-scale surveys of the matter distribution will be competitive with CMB 
constraints. However, in order to realise the potential of these surveys, we need to ensure that we are using a correct 
analysis, i.e. a general relativistic analysis, on cosmological scales. 

There are two fundamental issues underlying the GR analysis. 

• We need to correctly identify the galaxy overdensity A that is observed on the past light cone. The overdensity 
5 g defined in different gauges gives the same results on sub-Hubble scales, but leads to different results on large 
scales - and this remains true even if we use gauge-invariant definitions of S g . The observed A is necessarily 
gauge-invariant, and is unique. 

• We need to account for all the distortions arising from observing on the past light cone, including redshift 
distortions (with all GR effects included) and volume distortions. 

These GR effects come in to the measured 2-point correlation function in redshift space ^(rix, n 2 ; z±, z 2 )- We provide 
a new representation of this function which is relatively simple to calculate, and which takes into account all GR and 
wide-angle contributions. Our fully general relativistic wide-angle formalism recovers and generalizes previous work 
in the Newtonian plane-parallel (flat-sky) [TJ [2] and Newtonian wide-angle [SHSj cases. 



The GR analysis of the matter power spectrum [TUH2T| leads to corrections on cosmological scales of the standard 
Newtonian analysis (which is accurate on small scales). The observed galaxy overdensity is a function of the observed 
direction n and redshift z. It may be expressed in any chosen gauge. We use synchronous-comoving gauge, in which 



I. INTRODUCTION 



Observed galaxy density perturbation 




(1) 



In ACDM, we have W = [22"1 Eg] (a prime denotes d T ). 

We can write the observed overdensity at observed redshift z and in the unit direction n as 



A(n, z) = A s (n, z) + A K (n, z) + A 7 (n, z). 



(2) 
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Here A s is a local term (i.e. evaluated at the source) which includes the galaxy density perturbation, the redshift 
distortion and the change in volume entailed by the redshift perturbation. A K is the weak lensing convergence integral 
along the line of sight, and A/ is a time delay integral along the line sight. In the gauge ([!]), we have [T7] 
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Here xi z ) is the comoving distance, 5(z) is the bias and 

b e (z) = -(1 + z) 



dln[n g (l 
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(5) 
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where n g is the background number density. (We have changed some of the notation in [17 .) The directional 
derivatives are defined as 



d« = n J dj, d\ = tfdu, d?> = d\udi = d\\d\\ 



df l -2 X - 1 dn. 



(7) 

(8) 



The term <2(z) encodes the magnification bias, which arises from the perturbation to the flux of a galaxy, relative to 
a galaxy at the same observed redshift in the unperturbed universe (see |17)h 

The local term A s contains the Newtonian local terms, and in addition some GR corrections. The line of sight 
term A/ is a pure GR correction. The lensing term A K is the same as in the Newtonian analysis. Note that in 
we do not drop the radial derivative term — 2 X " 1 d^ which gives a negligible contribution to (W| on small scales, but 
should not be neglected on large scales. 

We can relate the metric perturbations to the matter density contrast in synchronous gauge, removing the residual 
gauge ambiguity. Taking into account that E" + aHE' — \^Gp m E — 0, we get 
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Here £l m (z) is the matter density and /(z) is the growth rate, 
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(9) 
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(11) 
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where D is the growing mode of 6. 



II. REDSHIFT-SPACE CORRELATION FUNCTIONS IN GENERAL RELATIVITY 

The correlation function can be decomposed in spherical harmonics as [13] 

£(ni,n 2 ,zi,z 2 ) = (A(ni, zi)A(n 2 , z 2 )) = C e (z 1 , z 2 )Y em (n 1 )Y e * m (n 2 ) . 



(14) 
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However, this turns out to be computationally expensive to implement for wide-angle correlations, and we follow 
the alternative decomposition used in previous analyses based on a Newtonian approach [BHS]- (The formula for the 
correlation function in redshift space, including wide-angle effects, was first derived in [3J[S], without GR corrections.) 
This alternative expands the redshift space correlation function using tripolar spherical harmonics, with the basis 
functions 
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(15) 



Here 



£i h h 
mi m 2 m 3 

is the Wigner 3j symbol. Because our triangles are closed, only two independent £'s appear. The expansion for £ 
may then be written in redshift space as a sum over these functions with L in the range \£i — £2] < L < l\ + 1% 
(see [6] for details). In the conventional expansion (14 1 of £ into CVs, the £ represents the angular momentum of the 
correlation function, and is summed to infinity. By contrast, in the tripolar expansion, the €s are representative of 
the radial derivatives du in the two directions ni and ri2. Consequently, the sums involved in the tripolar expansion 
are no longer infinite, making explicit computations in redshift space comparatively simple. 
The positions of two galaxies and their separation are given by (see Fig. [lj 



xi=xini, x 2 =X2n 2 , x i2 = xi - x 2 = Xi 2 ni 2 . 



(16) 




FIG. 1: Positions of a pair of galaxies on the lightcone. 
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The overdensity in the form ^ leads to the cross correlation functions in redshift space: 

£ab(xi,x 2 ) = (A^(xi)A B (x 2 )) = £ba(x 2 ,xi), A,B = s,k,I. 
First, we define the spherical transform of the matter overdensity [3] 

d 3 k 
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where Vg is a Legendre polynomial. Then by A s can be be decomposed as 
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(17) 
(18) 

(19) 

(20) 
(21) 

(22) 



Here a is a generalization of the Newtonian expression, /? has the same form as in the Newtonian analysis and 7 is a 
new term arising from GR corrections: see Fig. [2] below. 



In order to decompose £ ss , we need the correlator of (18 1: 

d 3 k 



(A^( Xl , Zl )A^(x 2 ,z 2 )) = (-1)* J ^( l fc)-(" 1 +" 2 )^ 1 (k.n 1 )^ 2 (k.n 2 )exp( i k.x 12 ) P s (k;z u z 2 ) 
where Ps(k; z\, z 2 ) is defined by 

(J(k 1 , 2l )<5(k 2 ,z 2 )) = (2^) 3 ^(k 1 +k 2 )P (5 (fc 1 ;z 1 ,z 2 ) . 
In terms of the primordial power spectrum and the transfer function, we have 

P s (k;z 1 ,z 2 ) = P prim (fc)T 2 (fc) J( ^g| Z2) . 

Expanding Ve and exp (ik • x) in spherical harmonics and applying the Gaunt integral [24], we obtain 
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This is an expansion in the tripolar basis functions (15), with coefficients 

£,2(x;zi,z 2 ) -- 



dk 
2^' 



^jUxk) Ps(k;z 1 ,z 2 ) 



Finally, we arrive at the tripolar decomposition of t; ss in the most general (GR wide-angle) case: 
^ s (x!,x 2 ) = b(z 1 )b(z 2 ) B ssn l2L (Xi,X2)Si 1 i 2L {n 1 ,n 2l n 12 )£,l{xi2]z ll z 2 ). 



(26) 



(27) 



(28) 



By ([19]), £ i: rii = 0, 1, 2 so that n = n\ + n 2 = 0, 1, 2, 3, 4; by ([26]) we also have L = 0,1, 2, 3, 4. The B coefficients for 
£ ss and the other £ab (see below) are given in Appendix [A] These coefficients involve the functions (20) and (22) that 
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contain GR corrections. In order to compare the GR and Newtonian cases, we set Q = for simplicity. We rewrite 
pOb as 
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(29) 



where N g is the comoving galaxy number density. Then for small redshift (i.e. x — > 0), we recover the Newtonian 
limit of a [2]: 
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In the same limit, (22 1 shows that 
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(31) 



which is nonzero. However, to recover the Newtonian limit, we have to take into account (27) when X12 0. Then 
we find: 



iCl -+ o, 

so that 7 drops out of £ ss in the Newtonian approximation. 



(32) 
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FIG. 2: The functions a,/3,y in (19 1, assuming a concordance model, and with Q = 0, b = 2 and n 9 as denned in [17] . 
unchanged by GR corrections, whereas 7 is not present in the Newtonian analysis. 



The remaining ah involve integrals along the lines of sight. The spherical transforms of A K , Aj are 



A K (n,z) = b(z) / dx<j(z,z) 



A 7 (n,z) = 6(z) / dx/i(«,«)^(x,l), 



(33) 
(34) 
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where 
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and x = x( z ): X = x(^)- These functions are illustrated in Fig. [3] a has the same form as in Newtonian analysis and 
\x corresponds to GR terms that vanish in the Newtonian limit. 
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FIG. 3: The functions a(3, z) and fi(3, z), assuming a concordance model, and with Q = 0, b — 2 and n g as defined in [17]. a 
is unchanged by GR corrections, whereas n is not present in the Newtonian analysis. 



For the lensing-lensing correlation, we find 

/Xl,X2 
dxidX2 B ^Kn t2L (Xi,Xi;X2,X2) Si 1 i 2L (n u n 2 ,n 12 )Q(xi2;z 1 ,Z2), (37) 

e lr l 2 ,L,n 

and for the 77 correlation 

/Xl,X2 
dx\dX2 - B //n l£2i (Xi,Xi;X2,X2)^ 1 ^L(ni,n 2 ,n 12 )^(xi2;5i,5 2 ). (38) 

The integration variables X\ 2 , &12 are given by 

Xi2ni2 = X12II12 + (xi - Xi) n i ~ (X2 ~ X2) n 2 , 

~2 -2,-2, XlX2 r 2 I ~2 , ~2M 
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XlX2 

Similarly, we find: 
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f K /(xi,x 2 ) = b(z 1 )b(z 2 ) / dxi^X2 Y • B K /n l£2L (xi,Xi;X2,X2)<S , £i£ 2 L(ni,n2,ni2)C£(xi2; z i,2 2 ), (43) 



where 
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(44) 
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Note that the remaining £ab follow from the symmetry in ( |17[ ). (See Appendix [A] for the explicit expressions.) 
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III. COMPUTING THE GENERAL RELATIVISTIC CORRELATIONS 

Here we consider the GR wide-angle correlation function £ ss (xi,X2) in various limits and then in the general case. 
For simplicity, and since it does not affect the comparisons, we assume Q = 0. 

GR small-angle (plane-parallel) limit 



For small angle and small galaxy separation, we have the plane-parallel or flat-sky limit, i.e. and n 2 are almost 
parallel [2 [2]. 

Using the properties of the Wigner coefficients and spherical functions, we get [5] 



Se 1 e 2 L{ni, ni, nn) 



tx h L 




Pi(ni-nxi), where ni ~ n 2 , nil = ni 2 



(46) 



With z\ ~ Z2 (xi ~ X2), Xn i s given by xn = xii n ii- From (28) we obtain the expansion of the redshift space 
correlation function in monopole, quadrupole and hexadecapole terms: 



£ ss (xi,xi) 
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(This is consistent with (TT1[T7].) In (47 1, we can divide the terms into 3 categories: those including only /3, those with 
a, and those with 7. The terms with p are the ones considered in the standard Kaiser analysis, and if we consider 
only them we will recover the classical results (see e.g. [5]). The terms including a are responsible for the "mode- 
coupling" effects, and they arise in the jacobian relating real- to redshift-space (for details, in Newtonian analysis, see 
e.g. [BHH])- These terms are usually ignored when using the plane-parallel approximation. Note that they vanish 
if the comoving radial distribution of galaxies is constant. Finally, the terms with 7 are the GR corrections. 

In Fig. |J we illustrate the GR corrections in ( |47[ ) to the Newtonian monopole and quadrupole of £ ss . For galaxy 
separations up to 600 Mpc, the GR corrections to the quadrupole remain negligible while the monopole is corrected 
by up to 2%. 
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FIG. 4: The GR corrections to the Newtonian plane-parallel limit of the £ ss monopole (left) and quadrupole (right), as a 
function of galaxy separation (Mpc). 



Note that in the Newtonian limit, a monopole and a quadrupole term proportional to a 2 survive: 



3x? 
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£ 2 2 7> 2 (ni -nn); 
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see also [3]. These terms are usually omitted from the standard flat-sky analyses [TJ [5] (for details see [5]). 

The remaining correlation functions which have not before been explicitly computed in the plane-parallel limit, 
are presented in Appendix |B| 



GR corrections on very large scales: analytical approximations 



For simplicity we consider the angular correlation, with z\ 
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where we impose a large-scale cutoff k mm , which we take as k min ~ Ho/2. We take X12 ~ Xh — 2(1 + z)H~ 1 (z). 
In this case, for k < l/xi2, we have P$(k;z) oc P pr i m (fc) = Ak ns , and jL{xi2k) — (xi2k) L / {2L + 1)!!. The second 
integral, for L > 0, can be approximated as 
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For L = 0, the integral vanishes because fc^ < 1/xn- Also, we can take P5 oc P pr im(&) since fct *C fc cq . Then, for the 
spectral index n s < 1, we obtain the analytical approximation: 
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(52) 



Note that, when n = 4 and L = 0, we require an infrared (IR) cutoff, fc m j n > 0, since £q becomes power-law divergent. 
(If n s = 1, there is a logarithmic divergence.) The IR cutoff appears only in the terms of the correlation function that 
contain Ylm with M = L — 0. (In this case Foo cx Vq = 1.) Therefore they add only an overall additive normalization 
to £ which is unobservable, and so the IR cutoff actually conveys no information. In addition, for large X12, the slope 
of is less steep when n increases. In other words, on large scales, B ss are bigger when they contain GR corrections. 
(The same is true for the other Bab-) 

Defining 9 as the angular separation of two galaxies, we have \12 = [2(1 — cos6')] 1 / 2 x- We assume that n g ~ Xii ' i 
where v is a suitable positive constant. Then by (201 and (22), in the large-scale limit x\2 3> Xh{z), we have 



a 
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Xh(z) 
b ( z )x 2 H (z] 
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(53) 
(54) 



Computation of GR wide-angle correlation function 

Finally, here we show the GR corrections to the full wide-angle correlation function. For illustrative purposes, we 
split the effects along the line of sight (n • k = 1), and transverse to the line of sight (n • k = 0). The line of sight 
case corresponds to a configuration where one galaxy is much further than the other and the angular separation 9 is 
small. In the transverse case, the two galaxies have the same redshift and 9 is not small. In Fig. [5] we plot the GR 
corrections for the two examples: 

• £ ss , as a function of z±, computed for z 2 — 0.1 and 9 = 0.1 rad. 

• £ ss , as a function of z\ = z%, computed for 9 = 0.3 rad. 

Along with the correlation function corrections, we plot the comoving galaxy separation, in order to have a clearer 
picture of the scales involved. In these computations we assumed a constant bias 6 = 2 and n g from the universal mass 
function prescription, see |17j . The plots of the GR corrections show the ratio of the correlation function computed 
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using the GR formalism with the one in the Newtonian case. Wide-angle effects are included in both cases, the only 
difference being the inclusion or not of terms deriving from 7. We decided to analyze two different cases, both for 
uniformity with previous literature and for physical reasons: pairs of galaxies (almost) along, and across the line 
of sight. The first case (left panel) corresponds to pairs that have a large difference in redshift and experience the 
maximum of the redshift-space distortions effect, while the second one (right panel) corresponds to pairs of galaxies 
with the same redshift, and for which the redshift space distortion effect is minimum. In a following companion paper 
we will present a detailed analysis of all the GR effects and their dependencies on scale, angular separation, bias and 
radial distribution of sources for all the terms presented in this work, along with predictions for observing these effects 
with future surveys. 
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FIG. 5: GR corrections to the Newtonian wide-angle £ ss correlation function (solid lines), along the line of sight (left; with 
Z2 = 0.1 and 8 = 0.1 rad) and transverse (right; with 8 = 0.3 rad). The corresponding comoving galaxy separation (Mpc) is 
also shown (dotted lines). 



IV. CONCLUSIONS 



We derived a fully general relativistic expression for correlation functions in the wide-angle case, that will need to 
be used in future surveys (such as SKA and Euclid) which measure galaxy correlations on very large scales. Our 
formalism recovers and generalizes previous results in the plane-parallel (flat-sky) and Newtonian approximations. 

We presented new results for the GR-corrected plane- parallel case. The GR corrections to the Newtonian monopole 
of £ SS) for galaxy separations up to 600 Mpc, are up to 2% (see Fig. [4]). 

Stronger GR corrections arise in the wide-angle case. We showed via illustrative examples that GR corrections on 
large scales (~ 1 Gpc separation) can in principle be of order 5-10%, as shown in Fig. [5j Of course, the observability 
of these effects is severely degraded by cosmic variance. Future large- volume surveys will reduce the cosmic variance - 
but the problem can in fact be removed if we are able to observe multiple tracers of the underlying matter distribution, 
as shown in |25j . In a companion paper we will present a detailed analysis of the predicted GR corrections for large- 
volume surveys, taking into account cosmic variance. 
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Appendix A: Coefficients in the tripolar decomposition of £ab 
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where a, = a(Zi), ft = and j t = 7(2*). 

For pfb: 
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K,K 1 



3V30 _~-l _ „ 
"5 X 2 ^ll<^22 > 



R 213 

D KK 1 



3V105 ~-l „ _ 
5 A2 ll a 22i 



121 _ 3V30 ~-l 



5 Xi ^li ^ 

= -3\/3 xr 1 ^ 1 °u CT 22 



r 123 _ 3V105 ~-l 

KK 1 _ 



R 112 

KK 2 



Xl O-ll^i 

3^ xr 1 ^ 1 0"il^25) 



where 0^7 = cr(zj, If 
For p8|): 



where = ju(2j, 
For fill: 



r 000 
°sl 2 



(A2) 



R 00° — 

n II 4 — Pll/ i 22' 



(A3) 



R 000 

-71M22, 



R 202 
fl s7 2 



^ /?1M22- 



(A4) 
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For (42): 



r> 000 

n SK 


= (1 + !j8i)*22. 




022 

o — 


VI (1 + o-aS. 


R Oil 

n SK 1 


= - 3>/3X2 1 (* + £&) "aS. 


B SK 


101 

1 — 


-\/3 xr 1 /3iai(r 2 2, 


R 121 

^SK 1 




B S« 


123 

1 — 




R HO 

fl SK 2 


= ^x^Vr^i^i^^^ 


B SK 


112 

2 — 




o 000 

fl SK 2 


= - 71 o- 2 2^ 


B SK 


022 

2 ^ 


-V5 71 ct 22 , 


R Oil 

n sK 3 


= 3a/3X2 1 7l°'22! 


B SK 


202 

o — 


-¥ /v 22) 


r> 220 

D SK 


= — "Jif /3lC 2 2, 


B sn 


222 

— 


21 A- ^) 


R 224 
fl ss 


_ 6N/70 o . 
— 35 Pl°22> 


B SK 


211 

1 — 


2V30~-lo _ _ 
5 X2 Pl cr 22' 



o 213 _ 2V105,~-1 

fl SK 1 — 



X 2 1 ^lO- 2 2- 



For (43): 



B kI 2°° — Cr ll/ X 22' 



R 202 
d kI 2 



-\/5 ^J^, 



Xl CTllM22- 



The remaining are 

/Xl 
dxi B i S n e2L (xi,Xi;X2)Se 1 e 2L (n 1 ,n 2 ,n i2 )(,l(xi 2 ;Si,z2), 

/Xi ^ 

dxi 2-j B K sn i2L (xi,Xi]X2)Si 1 i 2 L(n 1 ,n2,n i2 )£l(xi 2 ;Si 1 z 2 ), 

tl tl T „ 



,£ 2 ,L,n 



Xl,X2 



(A5) 



(A6) 

(A7) 
(A8) 



) = b( Zl )b(z 2 ) d Xl dx2 £ B /«« l£2i (Xi , Xi; X2> X2) 5^ 2 i(ni,n 2 , fii 2 ) CE(xi2; ^i, z 2 ), (A9) 



where 



Xi2 n i2 = (Xl - Xl) »i + Xi2ni2 

,2 _ ~2 1 ,,2 1 Xl r 2 /\,2 1 ,,2 



xf 2 = x? + xi + f [xl 2 -(xl + xl)], 



(A10) 
(All) 



and the corresponding B coefficients are 



^2° 00 = -(l +3^11, 



R 000 _ ... 
&Ia4 -72 Mil, 



B 



Is 3 

B Is 2 22 = ^2Mll) 



2«2X2 Mil; 



(A12) 



r 000 

n KS 



B KS \ ai 

R 211 
B ks 1 



r no 

D KS 2 



r 000 

n KS 2 



R 101 

"ks 3 



r 220 

D KS 



= (i + l&Ki, 

= 3^xr 1 ( 1 + 3^) CT n> 

= \/3 X^X^/^a^n, 
= - 72 on, 



r 202 
fl «s 



R Oil 

K KS 1 

R 213 

C KS 1 



R H2 



r 202 
d ks 2 



r 022 
d ks 



r 222 
D ns 



/5 (l + lA)^, 



^30 Xi *X 2 1 /?2a2CT li , 
-VE 72 (Til, 

2 V5 fl rr 
P 2 CT n , 



3 

2V70 
21 



R 224 

D KS 



-35- P20- n , 



r 121 _ 2y / 30 -,-1 o 
^1 - -5— Xl P20u, 



r 123 _ 2V105 ~-l a . 
-^s 1 5~ Xl P20-H • 



^Jk 2°° — Mll ^) 



r 022 
-°/k 2 



= -V5 



A*ll cr 22' 



R 011 

-O/k 3 



= 3V3 X 2 Vli CT 22- 



Appendix B: General relativistic £ab in the plane-parallel limit 



£ ss is given in (47). The remaining GR £ab have not before been given in the plane-parallel limit. 
Lensing-lensing correlation: 

/xi rxi ( 

dxv J dxv«riv<riv'\ 

3 f \^{xiv]z v ,z v ,) + — " — (,l{xvi")z v ,zy) J Vo{n x ■ n vl „) 

+ 3— — —tUxi'V, z v , z v ,)Vi(ni • ni'i") 



Xi'Xi'' 

+ *Ci&ixvi»\z v ,z v ,) — ' 2 

V 7 Xi'Xi" 



£ 2 (Xvi";zv,Zy) V 2 (ni ■ n vv >) 



\— — —^k{xvv';z v ,zyyP3{ixi ■ n V v>) + ^C4(xi'i";^i',^i")^4(ni • n V v>) 
Xi'Xi" 00 



II correlation: 



/Xi rxi 
dxv i dxi"fJ.ivfXii" $(xi'i";zi>,zy)'Po{n 1 -n V y). 

We assumed ni 2 = niiy. Then 

Xi'i"ni'i" = (xi' -Xi") n i' +Xi"Xr 1 Xnnn , 
Xi'i" = Xv +Xi" +Xi'Xv 

si correlation: 



/Xl 
dxi'Vu 



~ ( 1 + §^1 ) £a(XiV, z i, z v) + litoiXiv! z ii z v) 



-^■&{Xn';z u zi>)Vi(ni • n n .) + ^^(xiv, zi, z v )T 2 (ni ■ n lv ) 
Xi 3 
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where we assumed = nn/. Then 



Xivmv = (Xi -Xi') n i +Xi'Xi 1 Xunn 

2 2,2 

Xii' = Xi +Xi 



Is correlation: 



£ 7s (xi,xi) = b\ J dxiMii'j 



Xi 



1 + gft ) Co(Xi'i;^i',2i) +7i^ (Xi'i;^i',^i) 



(B6) 
(B7) 



"P (ni • ni/i) 



+ -^Ci(Xi'i;^i',^i)^i(ni • n ri ) + ^/?i£|(Xi'i;2i', 21)^2(111 • n vl ) 
Xi 3 



where we defined n;, = nyi. Then 



sk correlation: 



Cs«(xi,xi) = b l I dxv<?iv 



Xi'ini'i = -(xi -Xi') n i +Xii n n and xi'i = Xu' 



?y dXi'O-ii'l + $>{Xiv\zi,z v ) - +7i^o(Xii';^i^i') 



(B8) 



(B9) 



"P ( n i • nn/) 



+ 


K 


+ 




+ 


1* 



1 3 1 /3l«l \ ^1 , . 3 

— 1- i — pi + ^ Ciixn'^i^iO 7i?i(xii';^i>^i') 

:i' 5Xi' & xi / Xi' 

/3iai 



Xi 

1 

7 

[ — - — ) ^(Xii'^i^i')^^! • nn/) - ^^^(Xii'^i^iO^lni • nn/)} • 
VXi Xi'/ 35 J 



Xi'Xi 



7i ^(Xii'^i^i' 



"Pi(ni • nn-) 
"P 2 (ni • nn/) 



(BIO) 



ks correlation: 



/Xl 
dXi'O-ii 



1 + \&\ ) ^o(xi'i;^i',^i) - ( +71 ) £o(xi'i;zi',zi) 
5 / VXiXi' 



1 



3 I h Pi + \^°^ ) ^(Xi'i^iszi) - — 7id 3) (Xi'i;^i',^i) 

Xi' 5xi' 5 xi / Xi' 

1 - ^1) C 2 °(xi'i^i'^i) + (2-^- -71) d(xi'i;«i',«i) 
7 / V xi'Xi / 



"P (ni • ni/i) 
Pi(nj • ni/i) 
"P 2 (ni • ni'i) 



ai 2 
Xi Xi 



^^(Xi'i^i'^OPsCni-ni/i)- ^ft^(xi'i;*i',«i)P4(ni-ni/i)j . (Bll) 



Ik correlation: 



/Xi rXi 
dxv I dxi" juii'0-ii 



- Co(Xi'i";^i',2 ; i")^o(ni • ni/i») 



Xi' 



-^i(Xi'i»;^',2i«)Pi(ni • ni/i//) - (,1{xv\";zv,z vi )V2{thi • n-vv) 



(B12) 



kJ correlation: 



/Xi rXi 
c^Xi' / dxi" (TwVn 



Co 2) (Xi'i";2i',zi'')"Po(ni • n n //) 



H ^i 3) (Xi'i";-2i')^i")^i( n i • ni'i») -^ 2) (xi'i»;zi/,zi//)'P2(ni • n vv >) 

Xi 



(B13) 



Note that nyy> ^ —nyy xv = Xi"> an d n ii' 7^ — n i'i unless xi = Xi'- This is connected to the fact that we have 
already fixed galaxy 1 and galaxy 2 a priori and by the definition of nn, when we impose the plane-parallel limit. 
Indeed the second subscript 1 is the second galaxy and not the first one - i.e. in the flat-sky limit nn = — n 2 i. 
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